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При помощи аппарата неравновесной термодинамики сформулирован 
ячеечный метод рассмотрения явлений переноса в концентрированных дис- 
персных системах. В рамках этого метода равноправно учтены все неравно- 
весные процессы (вязкие, диффузионные и электрические} и получены вы- 
ражения для интегральных, т. е. наблюдаемых, термодинамических потоков 
и сил через локальные поля. Предложенный метод обеспечивает точное вы- 
полнение формулы Смолуховского для электроосмоса при использовании 
любой гидродинамической ячеечной модели (Кувабары, Хаппеля или 
Мехта). 


В отсутствие поляризации тонкого двойного электрического слоя 
(ДС) электрокинетические явления (электроосмос и потенциал течения) 
в дисперсии описываются формулами Смолуховского [1] 


9 а (1) 
Г [др—о 49 К 

СР (2) 
Ар [1-ю 4л9К 


где в, п, К — соответственно диэлектрическая проницаемость,. динами- 
ческая вязкость и электропроводность равновесного с дисперсией рас- 
твора, 5 — электрокинетический потенциал, $, — объемная скорость, 
Г — плотность тока, Аф и Ар — перепады потенциала и давления. Фор- 
мулы (1), (2) применимы к любой системе независимо от ее геометри- 
ческих характеристик (размера и формы частиц или капилляров, по- 
ристости) при 


ха» (3) 
Ве! — 2 (&*®_ Г]ма-<1 (4) 


где х — обратная дебаевская длина, а — размер частицы, фь = Еф ВТ — 
безразмерный штерновский потенциал. Неравенство (3) выражает усло- 
вие тонкости ДС, (4) — условие слабости его поляризации. 

Независимость формул (1), (2) от геометрии системы исчезает при 
переходе к поляризационному режиму Ке!->1, поскольку при этом на- 
рушается подобие электрического и гидродинамического полей в элект- 
роосмосе — основа вывода формул (1), (2). Поэтому. представляет ин- 
терес обобщение результатов теории поляризации тонкого ДС [1, 2], 
полученных для единичной частицы, на дисперсии с конечной объемной 
долей частиц. — 
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В подавляющем большинстве ситуаций, за исключением, пожалуй, 
только замечательного случая упомянутых выше двух электрокинети- 
ческих явлений в режиме Смолуховского, кинетические коэффициенты, 
описывающие процессы переноса, существенно зависят от геометрии си- 
стемы. Во всех этих случаях точное решение задачи о вычислении ки- 
нетических коэффициентов (КК), характеризующих концентрированную 
дисперсную систему, невозможно из-за очень сложной геометрии поро- 
вого пространства (даже когда форма дисперсных частиц предельно 
проста, например одинаковые сферы или цилиндры). 

Эффективным методом, в максимальной степени упрощающим зада- 
чу нахождения кинетических коэффициентов для концентрированной 
системы и в то же время учитывающим важнейшие особенности этой 
задачи, является переход к ячеечной модели. В этом методе рассмотре- 
ние концентрированной дисперсной системы, например, из одинаковых 
сфер, заменяется рассмотрением сферической ячейки с одной дисперс- 
ной частицей в центре. 

Размер ячейки определяется требованием, чтобы доля, которую со- 
ставляет объем частицы от общего объема ячейки, равнялась объемной 
доле дисперсных частиц в системе. Влияние остальных частиц эффек- 
тивно учитывается введением некоторых граничных условий на поверх- 
ности ячейки. 

При рассмотрении концентрированной системы на базе ячеечной мо- 
дели можно выделить три характерных масштаба: локальный, ячеечный 
‘и интегральный, которые соответствуют различным стадиям решения 
‘задачи. В локальном масштабе физически малый объем содержит мно- 
го молекул, но мал по сравнению с размерами частиц и с расстоянием 
между ними. В локальном масштабе мы задаемся кинетическими ко- 
эффициентами, характеризующими поровый раствор — такими, как вяз- 
кость жидкости и коэффициенты диффузии ионов. На этой основе мы 
записываем уравнения, описывающие распределения локальных термо- 
динамических сил и потоков в объеме ячейки, и соответствующие гра- 
ничные условия на поверхности частицы и на поверхности ячейки. За- 
тем возникает задача на базе определенных таким образом локальных 
распределений сформулировать выражения для «ячеечных» термодина- 
мических сил и потоков (т. е. таких, которые будут характеризовать 
ячейку как целое). Наконец, на последнем этапе мы должны выразить 
термодинамические силы и потоки, характеризующие систему в целом, 
через «ячеечные» характеристики и вычислить соответствующие кинети- 
ческие коэффициенты. 

Определение «ячеечных» термодинамических потоков и сил. Рацио- 
нальные выражения для термодинамических сил и потоков можно по- 
лучить [3], записав выражение для производства энтропии в системе 
(в нашем случае в ячейке). Производство энтропии в` ячейке получим, 
усреднив по ее объему выражение для локального производства энтро- 
пии, подобно тому как для другой системы это сделано в [4]. 

Запишем выражение для локального производства энтропии, следуя 
5] для случая отсутствия химических реакций, градиента температуры 
и объемной вязкости 


3 


: => до, : 
о + биз [| 6) 


Т Пе 9х, 9х; 


где Г — абсолютная температура, и. — химические потенциал моля ком- 
понента раствора &, #. — «диффузионный» поток компонента ©, Т. е. 
поток частиц относительно системы координат, связанных с каким-либо 
вотоком (в нашем случае —с потоком объема), Е. =—Р7„Уф — сила, 
действующая на частицы сорта , Пь — бесследовая часть тензора вяз- 
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ких напряжений (Пн==0), у — скорость жидкости, х: (1=1--3) — де- 
картовы координаты. | 

Выражение (5) определяет локальные феноменологические уравне- 
ния переноса: 


НИ в. о (6) 
А 
Фа = р № Це, — Ра) (7) 


ГДе 1, А« — локальные кинетические коэффициенты. 

При записи (6), (7) учитывался принцип Кюри ‚(потоки и силы раз- 
личной тензорной размерности не могут быть связаны между собой 
[5]), кроме того, в (7) отбрасывались члены, описывающие взаимодиф- 
фузию. 

Предположив идеальность раствора, (7) можно переписать в более 
удобном виде 


ЛоЮ 
У«= —— 


(у Са ++ = Сага, Уч ы —Да (ус + = сть УФ) (8) 


% 
где О. — коэффициент диффузии, 2. — валентность (с учетом знака 
заряда), с. — концентрация ионов, ф — электрический потенциал. 

Потоки (6), (7) удовлетворяют локальным уравнениям 


д 
`дх, (9 + Со:) =0 (9) 
| эк +0 (10) 
д; 
Е (11) 


где П’,=П,—Я’4, — полный тензор вязких напряжений, 9” — локаль- 


9ф 
Е Я Сага — объемная сила. 


-“ 

Уравнения (9) — (11) необходимо дополнить граничными условиями 
на поверхности, ограничивающей раствор. Для гидродинамически не- 
проницаемой и непроводящей частицы 

У|, =0 (12) 


(фм) |, =0 (13) 


где $ — поверхность частицы, п — нормаль к ней. Необходимо также 
задать условия на поверхности ячейки 51. Последние включают усло- 
вие, конкретизирующее гидродинамическую яЧеечиую модель: 

го у|„=0 (14) 


ное гидростатическое давление, |, = —Р 


(модель Кувабары [6]), 
(п. ПЬ) |.=0 (15) 
(модель Хаппеля [6], п,,  — нормаль и касательная к $:) 
9 059 |„=— 0; 51 6 |, (16) 
(модель Мехта и Морзе [7], где 9 — угол между нормалью п, и на- 
правлением макроскопического течения), и три граничных условия, 
отражающих условия измерения кинетических коэффициентов (на- 


пример, заданные значения трех потоков или трех сил). Чтобы кор- 
ректно сформулировать явный вид последних граничных условий, не- 
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обходимо выразить интегральные потоки и силы через ячеечные, а. их — 
через локальные. | 

Для решения этой задачи усредним производство энтропии (5) по 
объему ячейки 


орз) рыбы) 
| (17) 


При помощи уравнений (10), (11) и теоремы Гаусса второй член в (17) 
легко преобразуется в поверхностный интеграл 


И ‚ [Эра В о 
и 1. 38 (= — Ра) ола { шнаь | (18) 


$ 


Для удобства дальнейших вычислений разобъем потенциал ф на две 
составляющие 
ф=4Ф: + ф2 (19) 


выделив локально-равновесную часть $». Локальные концентрации 
ионов связаны с ф, уравнением Больцмана 


си — ®_ехр | —=") (20) 


где с(г) — локальная нормальная концентрация электролита. Выра- 
жение (20) является определением величин с(т) и ф». Возможность та- 
кой записи следует из основного постулата неравновесной термодина- 
мики — гипотезы о существовании локального равновесия [5]. 

Используя выражения (19), (20), уравнение (9) и формулу для 
химического потенциала ионов и„=ЮТ 1 с», объемный интеграл в (18) 
тоже можно преобразовать в поверхностный 


9: = — т |- Пу; Е ф1Р № 2а (фм -- бат!) + | У м 49: (21) 
$1 0. 0@=1,2 
где и= ЕТ шс(т); Пу = Пу — 96; Ф = 9'— ЕТ УМ — полное давле- 


0==1,2 


ние (с учетом осмотического). Локально-равновесный потенциал, как 
следует из его определения, не входит явно в выражение для произ- 
водства энтропии (21). Интегрирование ведется по внешней нормали 
к объему И; интеграл по поверхности $ равен нулю в силу граничных 
условий (12), (13). 

В отличие от работы [2], где рассматривается система пяти зави- 
симых потоков, будем рассматривать систему трех независимых пото- 
ков. Обозначим локальный диффузионный поток 


= УХ Ям (22) 
&—1,2 
локальный электрический ток 


й = у Его, [627 + соо») (23} 


@—1,3 


Введем «ячеечные» потоки: 
объема Г, 


‚В= = 9, 5ь (24) 
: 5% 
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‘электрический Г, 
А | { 45 (25) 
диффузионный Г 


те -- | | ре 49 (26) 


где $ — площадь экваториального сечения ячейки, 5$, — поверхность, 
натянутая на экватор ячейки. Тогда (21) можно переписать в виде 


3 
= —т У 1$» (27) 


где «ячеечные» градиенты давления Ф,, потенциала Ф, и химического 
потенциала Ф, определены в соответствии с '(21), (24) — (26) 


Ф, = — - ф Пия (28) 
У: 5 | 

— 1 

= — ; 29 
= $ фи (29) 
— 1 : 
Ф.; = 48 

з У, ТИ е (30) 


На основе (27) можно сформулировать «ячеечные» уравнения перено- 
‚са, которыми определяются «ячеечные» КК 9» 


= —У2„Ф» (31) 


Рть= ФР 


Введение интегральных потоков и сил. Между «ячеечными» и ин- 
ъ‘тегральными потоками нет принципиальной разницы. Определения 
«ячеечных» величин (24)— (26), (29) —(31) можно использовать и для 
определения интегральных величин в том случае, если получаемые 
«ячеечные» КК не противоречат двум следующим требованиям: 1) при 
выполнении условий (3), (4) и малости объемной доли « электропро- 
водность дисперсии должна выражаться формулой [8], 


Ф=К(1— Ка) (32) 


где ^=?/, для сфер, &=2 для ориентированных поперек поля цилинд- 
ров; 2) при выполнении условий (3), (4) отношения (1), (2) не дол- 
жны зависеть от объемной доли; в терминах уравнений (31) (1), (2) 
имеют вид: 


в _ 8 (33) 
58 АК 


Если же КК 5.» не удовлетворяют требованиям (39), (33), то мы 
введем новые определения интегральных потоков [», сил Ф, и уравне- 
‘ний переноса 


= — У ГиФь (34) 


т 


‘такие, что КК Г»„ будут удовлетворять требованиям (32), (33). Как 
будет показано в следующей работе, достаточно переопределить толь- 
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ко поток объема и соответственно интегральный градиент давления 


ЕЕ Ф 
1 = н@), Ф, =Н(@)Ф, (35) 
ь=1,  Ф,=Ф, (36) 
3= 1, Ф; = Фз (37) 


Такое переопределение обеспечит выполнение (32), (33) и одновремен- 
но сохранит вид выражения (27) для интегрального производства энт- 
‚ропии. Потоки и силы, найденные в соответствии с (35) — (37), явля- 
ются наблюдаемыми макроскопическими величинами. Чтобы однознач- 
но определить наблюдаемые КК из уравнений (34), необходимо зафик- 
сировать три из величин (35) — (37); тогда при помощи формул (24) — 
(26) и (28) —(30) можем получить три краевых условия на локальные 
поля, выполнение которых обеспечивает заданные значения зафикси- 
рованных величин. Эти три краевых условия совместно с условиями 
(12), (13) и одним из условий (14) —(16) позволяют найти однознач- 
ное решение уравнений (9)— (11), зная которое, можно по формулам 
(24) —(26), (28) —(30) вычислить значение трех искомых величин из 
фигурирующих в (35) — (37). 

Предложенная здесь процедура выражения наблюдаемых инте- 
гральных потоков и сил через локальные позволяет при исследовании 
электрокинетических явлений использовать любую из гидродинамиче- 
ских «ячеечных» моделей (14) —(16) и при этом получать непротиво- 
речивые результаты (в смысле выполнения структурной независимости 
формул (1), (2) при выполнении условий (3), (4)). Определенные в 
соответствии с этой процедурой интегральные КК удовлетворяют со- 
отношениям Онзагера, что позволяет для описания всех шести элект- 
рокинетических, диффузионно-электрических и диффузионно-кинетиче- 
ских явлений в концентрированных системах рещать задачу о нахож- 
дении только трех перекрестных кинетических коэффициентов. 
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ТНЕОКУ ОЕ МОМ-ЕООПЛВЕОМ ЕГЕСТЕВИСАТ ЗОВЕАСЕ РНЕМОМЕМА 
1№ СОМСЕМТВАТЕО РЕЗРЕВ$Е $У$ТЕМ$ 


1. АРРАСАТ!ОМ ОЕ ТНЕ МЕТНОР ОЕ ТНЕВМОБУМАМ!С$ ОР ВЕУЕКЗВЕЕ РАОСЕ$$Е$ 
ТО ТНЕ СЕМОГАВ МОРЕЕ ОР А СОМСЕКТКАТЕО Р1$РЕВ$1ОМ 


Звйоч У. М., Сваг ЕВ М. Г., ВогкочзЕаца Уи. В. 


Зит шагу 


Ву изше Че аррагафиз оЁ поп-едшИЬгиипй Шфегтодупапус$, опе паз 1огии!а{йеё Ше 
сей Шаг тео Гог сопз14егпя {Не рнепотепа о} фгапз{ег 1 сопсетшга{е4 @1зрегзе зуз{етз. 
\ИШт фе НатехуогК о! {515 шео4, опе Ваз едиаНу {фаКеп и соп$14егаНоп аП 1е поп- 
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{ог ИЦерта| (1. е., Фе оБзегуа Ме опез} Шегтодупаии!с Нихез ап@ Фогсез {тоивВ 1оса! 
Нез. ТНе зиасезей тео виагап{еез {Не ассигафе И тепё о! Не ЗтоиКНоузКу Гог- 
шша {юг е1ес{гоо$1031$ \Пеп изед Гог апу Вудгодупапис сеЙи]аг то4е| ({позе о? КиуаВа- 
га, Нарре!, ог Месв\). 


